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Our calculations have shown that the 6/5/2 — 7si/2 23131 cm~^ transition from the ground state 
in the ion is very sensitive to the temporal variation of the fine structure constant a = /he 



-75300 cm ). The line is very narrow, the ion has been trapped and laser cooled and the 



positive shifter line 5/5/2 — 6/7/2 4325 cm {q — -1-2900 cm ) may be used as a reference. A 
comparison may also be made with a positive shifter in another atom or ion. This makes Th^+ a 
good candidate to search for the a variation. 



PACS numbers: 31.30.Gs, 06.20.Jr, 31.15.am 



I. INTRODUCTION 

Theories unifying gravity with other interactions sug- 
gest temporal and spatial variation of the fundamental 
"constants" in expanding Universe (see e.g. review 
The spatial variation can explain fine tuning of the fun- 
damental constants which allows humans (and any life) 
to appear. We appeared in the area of the Universe 
where the values of the fundamental constants are con- 
sistent with our existence. The fundamental constants 
may be slightly different near massive bodies (see e.g. 
review 2] ) . There are some hints for the variation of dif- 
ferent fundamental constants in quasar absorption spec- 
tra [1, i, H, i, S B and Big Bang nucleosynthesis jl, 
data. However, a majority of publications report limits 
on the variations of the fundamental constants (see e.g. 
reviews [Hllil). 

The dependence of atomic transition frequencies on a 
may be presented in the following form 



w = wiab + qx, 
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In [1^ [T3| , it was proposed to use transitions with signifi- 
cantly different q factors for astrophysical and laboratory 
measurements of a variation. One can search for a varia- 
tion of a by comparing two frequencies of atomic transi- 
tions over a long period of time. Following the Ref. 
we can represent a measured quantity A{t) as 
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where Cj = doj/dt. Taking into account Eq. ([T]) we can 
rewrite Eq. ^ as follows 
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Narrow transitions with large and different q values are 
of experimental interest. 

Note that the atomic unit of energy cancels out in 
the ratio of two transition frequencies. The a depen- 
dence appears due to the relativistic corrections which 
rapidly increase with the nuclear charge Z, ~ Z^c?, 
and strongly depend on the electron angular momen- 
tum. Therefore, transitions in heavy atoms with larger 
electron angular momentum difference (like Al — 2 for 
s and d orbitals) have larger g-coefficients. At present 
the best laboratory constraint on the temporal variation 
of a of aja — (—1.6 ± 2.3) x yr"""^ was obtained 
by Rosenband tt al. in Ref. '16| by comparing the fre- 
quencies of the the ^£'5/2 transition in ^^^l{g+ 
{q = -52200 cm-i) and ^Sq ^Pq transition in 27A1+ 
{q — 146 cm~^). 

The relativistic corrections also rapidly increase with 
the effective charge ZcS which an external electron 
"sees". These corrections are proportional to Z^^. In 
the case of one electron above closed shells ZcS = Zi + 1 
where Zi is the ion charge. Therefore, highly charged ions 
are expected to have larger q. Unfortunately, the interval 
between the energy levels also increases ~ Z^g, therefore, 
there is a risk to be out of the laser range. However, the 
Coulomb degeneracy and configuration crossing phenom- 
ena may help here. Indeed, in a neutral atom an electron 
energy level with larger orbital angular momentum is sig- 
nificantly higher than a level with lower orbital angular 
momentum with the same principal quantum number n. 

For example, in the neutral Th the 5s electron is a 
core electron while the 5/ electron is a valence electron. 
Respectively, the one-electron energy of the 5/ electron 
is much higher than that of the 5s electron. Moreover, 
it is even higher than the energy of the 7s electron. As 
seen from the experimental spectrum of the energy levels 
of the neutral (four- valence) Th [l3|, the energy of the 
6d^7s5f state is higher than the energy of the 6d^7s^ 
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state. In the hydrogen-like Th the energy of the 5/ state 
is equal to the energy of the 5s state, i.e. it is significantly 
lower than the energy of the 7s state. Therefore, there 
should be an ion charge at which the level 5/ "crosses" 
the level 7s (at some higher charge it crosses 6s, etc). 
This kind of an approximate crossing (between the 5f6d 
and the 6d7s states) happens in Th^+ , where the interval 
between them is about 5500 cm~^ only. The interval 
between the 5/ and 7s states in Th'^+ is also relatively 
small. Finally, the lifetime of the excited state should be 
large to have a narrow line. 

All these requirements clearly point towards the 
5/5/2 — 7si/2 (23131 cm"-'^) transition in Th^+. It is also 
very important that the laser cooling of the 232rp]-^3+ j^j^ 
has recently been reported by Campbell et al. in their pa- 
per [3| ■ This was the first time when a multiply charged 
ion has been laser cooled. 

According to our calculations presented below the 
5/5/2 ~ 7si/2 transition is a negative shifter with q — 
—75300 cni^^. Another narrow line in the same ion is the 
positive shifter 5/5/2-5/7/2 (4325 cm'^) with q = -f 2900 
cm~^, this line may be used as a reference. Comparison 
may also be made with a positive shifter in another ele- 
ment, where q may exceed 30000 cm~^ (see the table of 
atomic clock transitions with q coefficients in review (l9| 
and recent work [20|). 



II. METHOD OF CALCULATION 

To find q factors we need to solve the atomic rela- 
tivistic eigenvalue problem for different values of a or, 
respectively, for different values of x from Eq. ([1]). The 
value of X was chosen to be equal to |a;| = 1/8. This is 
a convenient choice that allows us to neglect nonlinear 
corrections and, on the other hand, to make calculations 
numerically stable. Thus, we need to calculate atomic 
frequencies uj± for two values x ~ ±1/8. The corre- 
sponding q factor is given by 

g = 4(w+-tj_). (4) 

Since Th'^^ is a univalent ion we have carried out cal- 
culations of its energy levels in the frame of the Dirac- 
Hartree-Fock (DHF) method combined with many-body 
perturbation theory (MBPT). The latter allows us to 
take into account correlations between the valence elec- 
tron and the core electrons. 

We start from solving the DHF equations in the V^^^ 
approximation. On the first stage the electrons of the 
closed core were included in a self-consistency procedure 
and their orbitals were found. After that we constructed 
the valence orbitals for several low-lying states using the 
frozen-core DHF equations. The virtual orbitals were de- 
termined with the help of a recurrent procedure described 
in [2l|. The one-electron basis set of the following size 
was constructed: l-20s, 2-20p, 3-20d, 4-25/, 5-18^. 



TABLE L The low-lying energy levels (in cm~^) in the DHF 
and the DHF-I-E approximations are presented. The theoret- 
ical values are compared with the experimental data. 
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'The removal energy of the 5/5/2 state was found to be 
equal to 0.9414 au on the DHF stage and 1.0578 au on the 
(DHF-I-S) stage. The experimental value is 1.0588 au. 

At the next stage we included core- valence correlations 
(E) into consideration and the wave functions were de- 
termined by solving the equation 

H,«{E„)\^n)=E„\<Sj„), (5) 

with an effective Hamiltonian defined as 

H,ff{E) - Hfc + ^{E) , (6) 

where HpQ is the frozen-core DHF Hamiltonian and the 
self-energy operator E is the energy-dependent correc- 
tion, involving core excitations. In the following we will 
refer to it as the DHF+E formalism. 

III. DISCUSSION AND RESULTS 

In Table U we list the low-lying energy levels and com- 
pare them with the experimental data. To stress an im- 
portance of accounting for the core-valence correlations 
we present in Table [T] the results obtained on the stage 
of pure DHF approximation and in the frame of DHF+E 
formalism. 

As seen from Table |T] on the DHF stage even the or- 
der of the low-lying levels is incorrect. For instance, the 
6c?3/2 state lays deeper than the 5/5/2 state. An agree- 
ment between theoretical and experimental energy levels 
is rather poor. An accounting for the core- valence corre- 
lations (E corrections) recovers the correct order of the 
states. Besides that the theoretical energy levels become 
much closer to the experimental values. 

To find the q factors of the excited states in respect 
to the ground state 5/5/2 we need to carry out calcula- 
tions of frequencies lo± for two values x = ±1/8. These 
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TABLE II: The values of the q factors (in cm ^) found in the 
DHF+S approximations are presented. 
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calculations are similar to those carried out for the labo- 
ratory value of the fine structure constant aiab- For this 
reason we do not discuss them in detail. We only stress 
again that it is important to include the S corrections for 
obtaining the correct values of the frequencies lo± and, 
respectively, the correct values of the q factors. 

In Table im we present the q factors of the excited states 
listed in Table U in respect to the ground state, obtained 
on the DHF-fS stage. As follows from Eq. @ the ac- 
curacy of the q factors corresponds to the accuracy of 
calculations of uj^ and w_ . The spectrum of the energy 
levels of Th^+ is not too dense and we believe that the 
energy levels obtained at a-i- and a_ are found with the 
same accuracy as the energy levels computed at aiab- As 
seen from Table U the latter are reproduced within the 
10% accuracy. Correspondingly, the accuracy of the q 
factors can also be estimated at the level of 10%. 

As seen from Table HIl almost all q factors are negative. 
As expected, the largest q factor in absolute value was 
found for the 6/5/2 — 7si/2 transition. The only positive 



q factor was obtained for the 5/5/2 — 5/7/2 transition. 

The transition frequency between fine structure levels 
of one multiplct cijj_j_i (where J is the total angular 
momentum) in the first order in (aZ)'^ is given by the 
well known Lande rule: wj.j-i = AJ{aZY'. It directly 
leads to qj.j-i = (see Eq. (P)). 

A marked difference between 9(6/5/2 — 5/7/2) = 2900 
cm~^ and ^(6/5/2 — 6/7/2) = 4326 cm~^ demonstrates 
an importance of the second order relativistic corrections 
^ {aZ)'^ for Th. These corrections modify the expression 
for ujjj^i leading to (see, e.g., [22]) 

Luj.j-i = AJ{aZf -f {Bj - Bj^i) {aZf , (7) 

where A and Bj are certain coefficients and Bj are not 
small in comparison with A. 

Since Th is the heavy element with Z — 90, the param- 
eter (aZ) is not small for it. In particular, {aZ)'^ « 0.2 
and the term ~ ip-Zy^ gives a noticeable contribution to 
the 9(5/5/2 - 5/7/2). 

To conclude, we have calculated the q factors for 
a number of excited states in Th'^"'' in respect to the 
ground state 6/5/2- Our calculations showed that the 
9(5/5/2 — 7S1/2) is very large in absolute value and, re- 
spectively, the 5/5/2 — 7si/2 transition is very sensitive 
to the temporal variation of the fine structure constant. 
Since the 7s state is a metastable state, this transition 
is convenient for an experimental laboratory search of a 
variation. Another transition (6/5/2 — '^ji j-i) can be used 
as a reference. 
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